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Abstract 

Consider a generic triangle in the upper half of the complex plane 
with one side on the real line. This paper presents a tailored construc- 
tion of a discrete random walk whose continuum limit is a Brownian 
motion in the triangle, reflected instantaneously on the left and right 
sides with constant reflection angles. Starting from the top of the tri- 
angle, it is evident from the construction that the reflected Brownian 
motion lands with the uniform distribution on the base. Combined 
with conformal invariance and the locality property, this uniform exit 
distribution allows us to compute distribution functions characterizing 
the hull generated by the reflected Brownian motion. 

1 Introduction and overview 
1.1 Motivation 

Reflected Brownian motions in a wedge with constant reflection angles on 
the two sides were characterized by Varadhan and Williams in [EJ. Recent 
work by Lawler, Schramm and Werner [5] on SLE establishes a connection 
between one of these reflected Brownian motions (with reflection angles 
of 60° with respect to the boundary), chordal SLEg and the exploration 
process of critical percolation. They show that these three processes generate 
the same hull, using an argument that we shall outline in the following 
paragraph. The connection was carried even further by Julien Dubedat [IJ[3]. 
He compared SLEg and the aforementioned reflected Brownian motion in an 
equilateral triangle, started from a given corner and conditioned to cross the 
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triangle to a given point X on the opposite side. For these two processes, 
he showed that the conditional probability that the side to the right of the 
starting point is the last side visited before reaching X, is the same. He 
also showed that, under the assumption that SLEg is the scaling limit of the 
exploration process of critical percolation, this conditional probability can 
be used to prove Watts' formula for critical percolation (see [S]). 

The argument used in |2] to establish the connection between the three 
models is as follows. Let (Zt) be a stochastic process in an equilateral 
triangle T, started from a corner and stopped when it first hits the opposite 
side S. Define the hull Kt as the compact set of points in T disconnected 
from S by the trace of Zt up to time t. Denote by X the hitting point of S, 
and by r the hitting time. We shall refer to the probability distribution of 
the point X as the "exit distribution". For all three processes mentioned 
above, the exit distribution is uniform on the side of the triangle. Together 
with conformal invariance and the locality property (see Section 14.41 below 
and |S]), shared by all three processes, this exit distribution determines the 
law of the hull K T . The three processes therefore generate the same hull. 

This argument can be generalized to stochastic processes in arbitrary 
triangles that are conformally invariant and have the locality property: if 
the exit distributions of two such processes are the same, in particular if they 
are both uniform, then the processes generate the same hull. The papers of 
Lawler, Schramm, Werner ^| and Dubedat [3] show that in equilateral and 
isosceles triangles there exist reflected Brownian motions that have uniform 
exit distributions. The purpose of this paper is to generalize this result to 
arbitrary triangles, to review properties of these reflected Brownian motions 
and to discuss distribution functions associated with these processes and the 
hulls they generate. 

1.2 Notations and overview 

To present an overview of the present paper, we first need to introduce 
some notation. Given two angles a,/3& (0, it) such that a + (3 < it, we 
define the wedge W a p as the set {z G C : a — tt < argz < — (3}. We also 
define T a g as the triangle in the upper half of the complex plane such that 
one side coincides with the interval (0,1), and the interior angles at the 
corners and 1 are equal to a and f3, respectively. The third corner is at 
w a,p := (cos a sin /3 + i sin a sin /3)/ sin(a + /3). When a + (5 > 7T, the domain 
T Qj/ g is similarly defined as the (unbounded) polygon having one side equal 
to the interval (0, 1), and interior angles a and j3 at the corners and 1. We 
then identify the point at oo with the "third corner" w a ^. 
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Suppose now that and "&r are two angles of reflection on the left and 
right sides of the wedge W a $, respectively, measured from the boundary 
with small angles denoting reflection away from the origin ($L,f?R G (0, tt)). 
We shall use the abbreviation RBM$ L > $ R to denote the corresponding re- 
flected Brownian motion in the wedge W a 3. For a characterization and 
properties of these RBMs, see Varadhan and Williams ^3] (note that here 
we use a different convention for the reflection angles, namely, the angles i?i 
and $2 of Varadhan and Williams correspond in our notation to the angles 
■& L - tt/2 and $ R - tt/2). 

The main goal of this paper is to show that in every wedge W a 3 there 
is a unique RBM^ L $ with the following property: started from the origin, 
the first hitting point of the RBM of any horizontal line segment intersecting 
the wedge is uniformly distributed. This special behaviour is obtained by 
taking the reflection angles equal to the angles of the wedge, that is, = a 
and #r = (3. Restricting the wedge to a triangle we can reformulate this 
result as follows: 

Theorem 1 Let a, (3 G (0, tt), a + (3 < it, and let (Z t : t > 0) be an RBM a g 
in the triangle T a ^ started from w a ^ and stopped when it hits [0, 1] . Set 
t := inf{t > : Zt G [0, 1]} and X := Z T . Then X is uniform in [0, 1]. 

To prove this theorem we will cover the wedge W a p with a well-chosen 
lattice, and then define a random walk on this lattice. By construction, this 
random walk will have the desired property that it arrives on each horizontal 
row of vertices on the lattice with the uniform distribution. Taking the 
scaling limit then yields the desired result. The proof is split in two sections. 
In Section [2 we consider the easier case where the angles a and (3 are in the 
range (0, tt/2] such that tt/2 < a + (3 < tt. Section treats the extension to 
arbitrary triangles, which is considerably more involved. 

Section 21 collects some properties of the two-parameter family of RBMs. 
In particular, as was noted by Dubedat [1], we can show from the discrete 
approximations that the imaginary parts of the RBMs are essentially 3- 
dimensional Bessel processes. This allows us to describe the time-reversals 
of the RBMs, and sheds some light on how the RBM a 3 behaves in the 
domain T a ^ for angles a, (3 G (0, tt) such that a+f3 > tt. Finally, in Sectional 
we compute several distribution functions associated with the RBMs and the 
hulls they generate. This also reveals intriguing connections between RBMs 
started from different corners of the same triangle. 
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Figure 1: Picture of the lattice, showing the dimensions on the left, and the 
transition probabilities for a step of the random walk (from the origin in 
this picture) on the right. 

2 RBM in a restricted geometry 

Throughout this section, we assume that the angles a and (5 are fixed and 
restricted to the range (0, 7r/2] such that ir/2 < a + (3 < ir. We shall 
construct a random walk in the wedge W = W a ^ whose scaling limit is an 
RBM Qij a and whose horizontal coordinate is uniformly distributed all the 
time, and thereby prove Theorem ^for these restricted values of a and (3. 
The generalization to a generic triangle will be treated in Section El 

2.1 The lattice and further notations 

Throughout this paper we shall make use of a distorted triangular lattice, 
defined as follows. Let (p and ip be two angles in the range (0, 7r/2] such that 
7r/2 < ip + ip < 7r (as we shall see later on, the range for ip and ip is chosen 
such that the transition probabilities for our random walk are positive). We 
define as the set of vertices {j sin(<£> + ip) — k exp(i^) sin ip : j, k € Z} 
building the triangular lattice depicted in Figure ^ Throughout the paper 
we shall make use of the variables u := cosy sin ip, v := sin <p cos ip and h := 
sinysin^ to denote the lattice dimensions. When we use these variables, 
the values of ip and ip will always be clear from the context. 

To introduce some further notation, let us first consider how one defines 
a random walk (X n : n > 0) on that converges to standard complex 
Brownian motion in the full plane, before we consider the random walk in 
the wedge W in the following subsection. We set Xq := and for each 
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n > 0, X n is chosen among the nearest-neighbours of X n _i according to 
the probabilities a, b and c as depicted in Figure ^ We may then write the 
position X n of the random walk as a sum of steps S\ + S2 + • ■ ■ + S n where 
each step S n = U n + iV n is a complex- valued random variable taking on the 
possible values 

±(u + ih) with probability a; 
S n = ■[ ±(u + v) with probability b; (1) 
±(u — i/t) with probability c. 

To obtain a two-dimensional Brownian motion as the scaling limit of the 
random walk (X n ), it is sufficient that the covariance matrix of the real and 
imaginary parts U n and V n of each step is a multiple of the identity. This 
gives two equations for the probabilities a, b and c: 

(a + b) cot 2 ip + 2b cot 99 cot ^ + (b + c) cot 2 tj) = a + c, (2) 

a cot 99 — ccot V> = 0, (3) 

where cot x = 1/ tanx. The probabilities a, 6 and c can be determined from 
these equations, yielding 

a = A cot ^(cot 99 + cot ip) , (4) 
6 = A(l — cot tp cot ip), (5) 
c = A cot ip (cot 99 + cot ip), (6) 

where A = ^[cot </?(cot 99 + cot tp) + sin -2 is the normalization constant. 
One may verify that ip and ip must satisfy 7r/2<p + ilj<iTto make all 
three probabilities nonnegative. 

We conclude this subsection with a short discussion of how one obtains 
the scaling limit of the random walk (X n ). To do so, for every natural num- 
ber N > one may define the continuous-time, complex-valued stochastic 
process {z[ N ^ : t > 0) as the linear interpolation of the process 

Y t (N) = j^X lN . tl (7) 

making jumps at the times {k/N 2 : k = 1,2,...}. Here, a 2 is the variance of 
the real and imaginary parts of the steps S n , that is, a 2 = E[U 2 } = E[V 2 ]. 
It is then standard that z[ N ^ converges weakly to a complex Brownian 
motion in the full complex plane C when N — > 00 (topological aspects are 
as described in the Introduction of Varadhan and Williams |15j). 
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Figure 2: Definition of the reflected random walk in a wedge. 



2.2 Reflected random walk in a wedge 

We now return to the wedge W = W a p with angles a and (3 fixed in the 
range (0, 7r/2] such that ir/2 < a + j3 < it. Clearly, this wedge is covered 
nicely with vertices of the lattice when we set cp := a and ip := /3, see 
Figure |2j For the duration of this section we consider these values of ip and tp 
to be fixed. Later, when we generalize to arbitrary triangles, the relation 
between ip, ijj and a, f3 will not be so simple, which is why we already reserve 
the symbols cp, if) to denote the angles of the lattice. 

We shall denote by G = G a $ the set of vertices obtained by taking the 
intersection of T = T^^ with W. We shall call the vertices of G having six 
nearest neighbours along the lattice directions interior vertices. The origin 
will be called the apex of G, and the remaining vertices will be referred 
to as the boundary vertices. The set of boundary vertices may be further 
subdivided into left boundary vertices and right boundary vertices, with the 
obvious interpretation. 

Given a vertex x of G, a reflected random walk {X^ : n > 0) on G 
is defined as follows. We set X§ := x and for each n > 0, if X*_i is an 
interior vertex, then X£ is chosen among the six nearest-neighbours of X^_ 1 
according to the probabilities a, b and c as before. This guarantees that the 
scaling limit of the random walk is Brownian motion in the interior of the 
wedge. It remains to specify the transition probabilities for the random walk 
from the boundary vertices and the apex. 

In this paper, we restrict ourselves to the case where the transition prob- 
abilities on all left boundary vertices are the same, and likewise for the right 
boundary vertices. We write E 2 [5i] for the expected value of the first step 
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of the random walk started from the vertex z, and assume that it is nonzero 
at all the boundary vertices and the apex. Then, as we shall prove in Sec- 
tion the random walk converges to a reflected Brownian motion. More- 
over, the directions of reflection on the sides of the wedge are given by the 
directions of E z [Si] at the left and right boundary vertices. It follows that 
by playing with the transition probabilities from the boundary vertices we 
can obtain different RBMs in the scaling limit. 

There is, however, only one choice of transition probabilities for which 
the random walk has the following special property: started from the origin, 
the random walk first arrives on any row of the lattice with the uniform 
distribution on the vertices of that row. We will refer to this special case 
as the uniform (random) walk. To derive its transition probabilities one 
proceeds as follows. In a picture where we represent the steps of the walk 
by arrows, we have to make sure that every vertex in a given row has two 
incoming arrows with probability b from vertices in the same row, and two 
incoming arrows with probabilities a and c from vertices in the rows above 
and below. This completely determines the transition probabilities from the 
boundary vertices, see Figure El for a picture of the solution. 

In formula, if x is a left boundary vertex, then we have the following 
transition probabilities for the uniform walk: 



and the transition probabilities from the right boundary vertices are defined 
symmetrically, as shown in Figure EJ At the apex we simply choose the 
transition probabilities to each the two vertices directly below the apex 
equal to 1/2. 

It will now be convenient to decompose the position X® at each step 
of the uniform walk as J n (u + v) — K n (u + ih). Then K n is a nonnegative 
integer denoting a row of vertices on the lattice, and J n is a nonnegative 
integer denoting the position on the K n ih row. Observe that there are a 
total of N(k) = k + 1 vertices on the kth row, so that J n ranges from to 
N{K n ) — 1. Henceforth, we shall always adopt this convention of numbering 
rows on the lattice in top-down order, and vertices on each row from left to 
right. 

By construction, the uniform random walk (X®) started from the origin 
has the property that at every time n > 0, the position of the walker is 



p[x, x] = p[x, x + (u + v )] = b, 
p[x, x + u + ih] = a, 
p[x, x — u — ih] = a + c, 
p[x, x + v — ih] = c, 



(8) 
(9) 
(10) 

(11) 
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uniformly distributed on the rows of the lattice. More precisely, for the 
uniform walk (X®) the following lemma holds: 

Lemma 2 For all n E N, if ko,k\, . . . ,k n is a sequence of natural numbers 
such that ko = and \k m — k m -i\ < 1 for all m = 1,2, ... ,n, then for each 
j = 0,l,...,N(k n )-l, 

P[J n = j | Kq = ko, . . . , K n = k n ] = P[J n = j I K n = k n ] 



N{k n 



In particular, if k E N and T := min{n > : K ri G k} is the first time at 
which the walk visits row k, then for each j = 0, 1, . . . , N(k) — 1, 



N(k) 

Proof: The first claim of the lemma is proved by induction. For n = 0, 1 
the claim is trivial. For n > 1, 

■P[«^7i = J I -^n = ^n] 

S P[J"n = i,-^n = K | #n-l = K-l] P[-^n-l = fcn-l] 

= hzl (\2\ 

E P[-K"n = k n I K n _i = k„,-i] ¥[K n -i = kn-l] ' ^ ' 

Using the induction hypothesis, it is easy to compute the conditional prob- 
abilities in the numerator and denominator of this expression from the 
transition probabilities given earlier. We first assume k n > 1. Then for 
k n -l = k n these conditional probabilities are 2b/N(k n ) and 26, respectively. 
For k n -i = k n ± 1 they are (a + c)/N(k n -i) and (a + c)N(k n )/N(k n _i), 
and for all other k n -\ they vanish. When k n < 1 the contributions from the 
apex take a different form, but the computation is equally straightforward. 
Now observe that the quotient in Equation (j!2j) has the same value if the 
event {K n -\ = A; n _i} is replaced by {Kq = ko,...,K n -i = k n -\}. The 
first claim of the lemma follows. The second claim of the lemma is then a 
straightforward consequence. □ 

To prove convergence of the uniform walk to a reflected Brownian mo- 
tion, we shall need an estimate on the local time spent by the walk on the 
boundary vertices of G. This estimate is provided by Lemma |31 below. The 
proof of Lemma 01 is postponed to Section 14.11 since the proof will give us 
some intermediate results that we will need only in Section |1J 
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Lemma 3 Consider the walk (X®), killed when it reaches row M > 0, and 
let x be a fixed boundary vertex of G. Then the expected number of visits 
of (X®) to x before it is killed is smaller than (a + c) _1 . 



2.3 Identification of the scaling limit 

Given a point x in W, the scaling limit of the reflected random walk started 
from a vertex near x is obtained in a similar way as described in Section \2. II 
That is, if (xn) is a sequence in G such that \xn — Nax\ < 1/2 for all 
N € N, then for every natural number N > we define the continuous 
complex-valued stochastic process (z[ N ^ : t > 0) as the linear interpolation 
of the process 

Y > m = i^^j < 13 > 

making jumps at the times {k/N 2 : k = 1,2,...}. Here, a 2 is again the 
variance of the real and imaginary parts of the steps S n = X n — X n -\ of 
the random walk in the interior of the wedge. Then it is clear that the 
random walk converges to a complex Brownian motion in the interior, but 
the behaviour on the boundary is non-trivial. 

It is explained in Dubedat [I] how one proves that the scaling limit of 
the random walk (X®) is in fact reflected Brownian motion with constant 
reflection angles on the two sides, started from x. The core of the argument 
identifies reflected Brownian motion as the only possible weak limit (Zf) 
of (Z^). We repeat this part of the argument below, mainly to allow us 
in Section 01 to point out some technical issues that need to be taken care of 
when we generalize to arbitrary triangles. 

We use the submartingale characterization of reflected Brownian motion 
in W (see |15j . Theorem 2.1), which states the following. Let C 2 (W) be the 
set of bounded continuous real- valued functions on W that are twice contin- 
uously differentiable with bounded derivatives. Assume that the two reflec- 
tion angles #£ and #r on the sides of W are given. Then the RBM^ L) ^ H in W 
started from x 6 W is the unique continuous strong Markov process (Zf ) 
in W started from x such that for any / € C 2 (W) with nonnegative deriva- 
tives on the boundary in the directions of reflection, the process 
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f(Z*)-- I Af(Z*)ds, (14) 



o 



where A is the Laplace operator, is a submartingale. 

To prove convergence of our random walk to an RBM, it is therefore 
sufficient to show that there are two angles "&i and "&r such that if / is a 
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function as described above, then for all < s < t 



liminfE^ 

N^oo 



f(z t ) - f(z s 



Af(Z u )du 



> 0. 



(15) 



Here, Ej^ denotes the expectation operator for the Nth approximate pro- 
cess started near x, and we have dropped the (N) superscript on (Z t ) to sim- 
plify the notation. Note that the conditioning in the submartingale property 
is taken care of because the starting point x is arbitrary. In fact, it is suffi- 
cient to verify (|15j) up to the stopping time t := inf{i > : ImZt < — M} 
for some large number M. We will make use of this to have uniform bounds 
on the error terms in our discrete approximation. 

We now show (|15(l . Let / be a function in C 2 (W), and write u k = k/N 2 
for the jump times of (Zt). Then by Taylor's theorem, there exist (random) 
times T k between u k and u k+ \ such that 

f(Z Uk+1 ) - f(Z Uk ) = -±-(d x f{Z Uk )U k+l + d y f{Z Uk )V k+l ) + 

^±-^(d 2 f(Z Tk )U 2 +1 + d 2 f(Z Tk )V 2 +1 + 2d x d y f(Z Tk )U k+1 V k+1 ), (16) 

where U k and V k are the real and imaginary parts of the fcth step of the 
underlying random walk, as before. 

When we now apply the expectation operator, we distinguish between 
the behaviour on the boundary and in the interior (we may ignore what 
happens at the apex because of the negligible time spent there, see below). 
If B denotes the set of boundary vertices of G, then we may write 



EM [f(Z 



Uk + 1 - 



f(Z Uk )]=Bl 



(N) 



2N 2 



A/(Z r .) + 



(j^d x f(Z Uk )U k+1 + ±-d y f{Z Uk )V k+l + 0(N- 2 ] 



^-{NaZ u .eB} 



(17) 



where the 0(N 2 ) error term is uniform in {z E W : Im z > — M }. Summing 
over k from [N 2 (s A r)J to [N 2 (t A r)J yields 

Ef)[/(Z tAT )-/(Z iAT )]=Ef) / -A/(Z„)d« +o(l) + 

S E ^ {iv^ /( i\fe )E ^ l] + lt dyf{ lk )E > [Vl] +°^ 2 )} -( 18 ) 

z<=B 

In the last line, L z is the local time at z (the number of jumps of the random 
walk to z) between times s At and t At, and E z is expectation with respect 
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Figure 3: The thick arrows in this figure represent the directions of reflection 
of the reflected Brownian motion described in the text. The reflection angles 
are such that the Brownian motion will hit the dashed line with uniform 
distribution. 



to the underlying random walk started from z. Lemma 01 shows that when 
the starting point of the walk is 0, the expected local time spent on the 
boundary vertices up to time r is of order O(N). Then one can use the 
Markov property to see that this is enough for us to ignore the 0(N~ 2 ) 
term in the limit. 

Observe that the remaining term in braces is just the derivative of / 
at the boundary point z/Na along the direction of E 2 [Si]. Thus, if the 
function / has nonnegative derivatives along this direction on the boundary 
of W, then the term in braces in equation (|18|) is nonnegative, and the 
desired result (jl5fl is obtained. This proves that the scaling limit of the 
reflected random walk is a reflected Brownian motion, and that the angles 
of reflection with respect to the left and right sides are given by 

^ L + a = arg(-E z [5 1 ]) and R + p = - arg (E y [5i]) , (19) 

where z is any left boundary vertex and y any right boundary vertex. It is 
clear from the computation that this result is valid generally for any choice 
of transition probabilities on the left and right sides. We shall however focus 
again on the special case of the uniform walk introduced in Section 12.21 

The computation of the reflection angles for the uniform walk is straight- 
forward. By symmetry it is enough to compute only From the transition 
probabilities one may verify that 

. f a , fl x E z [ReSi] (q - b) cot ip - (a + b) cot ip 
cot(^ + P) = -^^J = Ya • (20) 



Substituting equations then yields 



cotO? fl + /3) = "ft . 1 = cot(2V) = cot(2/3). (21) 
2 cot w 



11 



Thus, the angle of reflection of the Brownian motion with respect to the 
right side is simply $r = [3. Similarly, the angle of reflection on the left side 
is given by = a. 

We conclude that the scaling limit of the uniform random walk defined 
in Section 12.21 is a reflected Brownian motion with fixed reflection angles 
on the two sides of the wedge. The angles of reflection are a and (3 with 
respect to the left and right sides, respectively, as illustrated in Figure El 
It follows from Lemma that the RBM a ^ has the special property that in 
the wedge W a ^, the RBM first arrives on any horizontal cross-section of the 
wedge with the uniform distribution. By a simple translation it follows that 
in the triangle T a ^, the KBM a ^ started from the top w a ^ will land on [0, 1] 
with the uniform distribution. This completes the proof of Theorem ^ f° r 
angles a, (3 € (0,7r/2] that satisfy tt/2 < a + (3 < tt. The case of arbitrary 
angles will be considered in the following section. 

3 RBM in a generic triangle 

In the previous section we proved Theorem ^ for angles a, (3 € (0, 7r/2] such 
that 7r/2 < a + (3 < tt. Here we generalize to generic triangles. That is, 
we now choose the angles a and f3 arbitrarily in the range (0, tt) such that 
a + (3 < tt. These values of a and (3 are assumed fixed for the remainder of 
this section. As in Section|21 we will define a uniform walk in the wedge W at /3 
and identify the scaling limit as a reflected Brownian motion. 

3.1 Choice of the lattice 

The first thing that we have to do is to find a lattice covering the wedge 
W = W a: p in a nice way. We will make use of the triangular lattices ^ 
introduced in the previous section. Consider such a lattice T^^, and choose 
two vertices xl and xr on the first row of the lattice such that xr is to the 
right of xl. Then the two half-lines {txi ■ t > 0} and {txj? : t > 0} define 
the left and right sides of a wedge that is covered nicely by vertices of ^ . 
We will show below that for any given wedge W there is a choice of the 
lattice angles (p,ip and the two points xl and xr such that the wedge thus 
defined coincides with W. 

We start by introducing some notation. Given the lattice T v ^ and the 
two vertices xl and xr, we may introduce two integers ni and ur that count 
the positions of xl and xr on the first row of the lattice. More precisely, 
the integers til and ur are defined such that xl = —u — ni{u + v) — ih and 
xr = v + tir{u + v) — ih (for an illustration, see FigureEJ). Conversely, given 
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n R = -2 
or til = 1 



tir = -1 
or til = 



n R = 
or rii = — 1 



Figure 4: Different wedges can be covered by the same triangular lattice, 
by changing the directions of the two sides as shown. These directions can 
be expressed in terms of integers til, and ur, as explained in the text. Each 
thick line in the figure can represent either the right side of a wedge (for 
which the corresponding value of ur is given in the figure), or the left side 
of a wedge (for which the corresponding value of ni is given). 

two integers ni and ur and the lattice the vertices xl and xr are 

fixed by these equations. Observe that til and ur must satisfy ni + tir > 
to make sure that xr lies to the right of xl- Figure 21 shows the wedges one 
obtains for different choices of til and ur on a given lattice. 

The main claim of this subsection is that for any choice of the angles a 
and (3, there is a choice of integers hl,ur and of the lattice angles (p,ip 
such that the wedge one obtains as described above coincides with the 
wedge W a $. This results is a direct consequence of Lemma @] below. The 
proof of the lemma, which will give explicit formulas for til, Ur and <p, ip in 
terms of the angles a and (3, is postponed to the end of this subsection. 

Lemma 4 Let a, (3 G (0, it) such that a + (3 < tt. Then there is a choice 
of (possibly negative) integers til and ur with + ur > 0, and angles 
(fjifi G (0, 7r/2] with 7r/2 < ip + ip < tt, such that 



We remind the reader that in Section |21 we considered a random walk on 
a wedge with angles a, j3 G (0, ir/2] such that ir/2 < a+(3 < tt. The situation 



cot a = ni(cot (p + cot ip) + cot ip; 
cot (3 = nft(cot ip + cot ip) + cot ip. 



(22) 
(23) 
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described there corresponds to the special choice of til = n R = 0, <p = a 
and ip = /3 in Lemma HJ Thus, the random walk construction of Section [21 
is a special case of the more general construction we are considering in this 
section. We would also like to remark at this point that the choice of ni, n>R 
and tp, tp in Lemma 0] is not unique in general. For instance, we will see in 
the proof of Lemma 21 that it is always possible to choose the angles tp, ip in 
the range [7r/4, n/2]. Thus, there are angles a, (3 such that the choices made 
in the proof of the lemma and in the construction of Section |2] are different. 

From now on, we will assume that the values of ip and ip are fixed as 
in Lemma 0J so that the lattice T = is fixed. As we explained above, 
the lattice provides a nice covering of the wedge W = W a a- Following the 
conventions introduced in Section EJ we shall denote by G = G a $ the set of 
vertices obtained by taking the intersection of T with W. We again subdivide 
the set G into the apex, interior vertices and left and right boundary vertices. 
We now conclude this subsection with the proof of Lemma 0] 

Proof of Lemma IH Assume first that both a and (3 are smaller than \. 
Then we can take 

til = [cot a] — 1, tir = [cot/3] — 1, (24) 

and solve equations (|2*2*)) and ((2*5|) for ip and ip to obtain 

n R + 1 n L 

cot ip = cot a cot p; (25) 

n L + n R + l n L + n R + l 

cot ip = cot (3 — cot a. (26) 

nL + ri R + l n L +n R + l 

Observe that since ni < cot a < ul + 1 and ur < cot j3 < tir + I, the angles 
ip and ip are in the range [ir/A,ir/2]. 

It remains to consider the case when either a or (3 is at least 7r/2, and 
by symmetry it suffices to assume a > ir/2. Then, we can for instance 
set k = [cot q + cot /3] , and let I be the smallest positive integer such that 
l/(k + I) > — cot q/ cot /3. We then set 

n L :=-l, n R :=k + l-l, (27) 

and the angles ip and tp are given by the equations (|23|) and (|2li|) as before. 
From the inequalities 

I cot a I — 1 , 
> — > r—r 7 , (28) 



k + l cot /3 _ k + I - 1 
plus the fact that k > cot a + cot /?, it follows that tp and ip are in the 
range [7r/4,7r/2]. This completes the proof. □ 
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3.2 Discussion of the random walk construction 

Now that we have chosen the lattice on the wedge W, the next task is to 
define a random walk (X*) on G whose scaling limit is reflected Brownian 
motion. As in Section [2.21 we will focus on the special case of the uniform 
walk, i.e., the random walk on the lattice that stays uniform on the rows 
all the time. Earlier we explained how one defines this random walk in the 
case til = ur = 0. In this subsection we will describe what one has to do 
to generalize to other wedges, i.e., to the case where ni or hr or both are 
nonzero. In the following subsection we will then spell out the transition 
probabilities of the uniform walk for this general case. 

It is clear that we should define the transition probabilities from the inte- 
rior vertices in the same way as before. This will guarantee that the random 
walk will converge to Brownian motion in the interior of the wedge. The 
nontrivial task is to define the transition probabilities from the boundary 
vertices and at the top of the wedge. Once these are defined, we will use the 
strategy of Section 12.31 to identify the scaling limit as a reflected Brownian 
motion. Let us therefore reconsider the arguments used there to identify the 
scaling limit, to see if they still apply to general wedges. 

Looking back at Section 12.31 we see that it is still sufficient to verify 
Equation Q15|) for appropriate functions /. Moreover, Equation Q18JI is still 
valid after making a Taylor expansion. However, we emphasize that in the 
previous section the sum over the boundary vertices picked up contributions 
from only one left boundary vertex and one right boundary vertex on every 
row of the lattice. This is no longer the case in the general situation, as we 
shall now discuss. 

Remember that the boundary vertices are defined as those vertices hav- 
ing less than six nearest-neighbours in G. Now consult Figure |1J Then 
we see that for a given value of ul, the number of left boundary vertices 
on each row of the lattice is fixed and equals Nl = \ul\ + ^{n L >o}- Like- 
wise, for given ur the number of right boundary vertices on every row is 
Nr = \ur\ + l{ nil >o}. (An observant reader will notice that if til or ur is 
negative, then on the first few rows of the lattice the number of boundary 
vertices is less than Nl + Nr. This is a complication that we will deal with 
later on.) 

We conclude that the sum in Equation (JT8|) picks up Nl contributions 
from the left boundary vertices in every row. Each of these contributions 
is proportional to the derivative of / near the boundary in the direction of 
E z [Si] . We have to guarantee that these contributions add up to something 
positive if / has positive derivative along some (yet to be determined) fixed 
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direction on the left side of the wedge. To do so, it is clearly sufficient to 
impose that E^fSi] is the same at all the left boundary vertices, and this is 
what we shall do. The corresponding condition will be imposed at the right 
boundary vertices. 

It is clear from the discussion in Section 12.31 that the scaling limit of the 
random walk will be a reflected Brownian motion. Moreover, the directions 
of reflection on the two sides of W will be given by the values of E 2 [Si] at 
the left and right boundary vertices. Let us remark that the condition that 
E z [Si] is the same at all the boundary vertices may be more than necessary, 
but by imposing this condition we avoid having to estimate the differences 
in expected local time at different boundary vertices. 

In summary, we will look for transition probabilities from the boundary 
vertices that satisfy the following conditions: 

1. The summed transition probability to a given vertex from vertices in 
the row above is a + c, and likewise from vertices in the row below. 

2. The summed transition probability to a given vertex from vertices in 
the same row is 26. 

3. The expected value of the first step of the walk from every left bound- 
ary vertex is the same, and likewise for the right boundary vertices. 

Observe that condition 1 introduces an up-down symmetry which is not 
inherent in the geometry of the problem, but will be of importance later. 
Explicit expressions for all the transition probabilities of the uniform walk 
will be given in the following subsection. This will show that for any ni, 
it is possible to choose the transition probabilities such that they satisfy 
conditions 1-3. 

To conclude this subsection, we point out one further complication that 
we already commented on earlier. This complication arises when either nx 
or riR is negative, because then the first few rows of the lattice will contain 
less than Nl + Nr vertices in total. Hence we will have to specify separately 
the transition probabilities on these first rows of the lattice. As we shall 
discuss below, this complication is quite easy to handle. 

Consider once again Figure 0] and recall the definition of the integers ni 
and hr in Section 13.11 illustrated in the figure. Then one notes that each 
row of the lattice contains exactly til + npt + 1 vertices more than the row 
above. In other words, the total number of vertices on row k of the lattice 
is N(k) = {til + rifl + l)k + 1. From this one can compute that if either hl 
or tir is negative, then the first row of the lattice that contains at least all 
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of the Nl + Nr boundary vertices is row ko, where 

\ul - n R \ 



n L + n R +l 



(29) 



For example, for ni = 2 and tir = — 1 the numbers of left and right boundary 
vertices are Ni = 3 and Nr = 1, respectively. However, the first two rows 
of G contain only -/V(O) = 1 and -/V(l) = 3 vertices in total, so that the first 
row of the lattice that contains at least Nl + Nr vertices is row 2. 

The remaining problem is then to define the transition probabilities from 
the vertices in the first ko rows of the lattice. It is not very difficult to choose 
the transition probabilities on these first few rows such that the walk will 
trivially stay uniform on the rows of the lattice. In fact, we can choose the 
transition probabilities such that conditions 1 and 2 above are also satisfied 
at rows 1 up to ko — 1, as we shall see below. This guarantees that the walk 
does indeed stay uniform on rows. We can not make condition 3 hold on the 
first ko rows. This is no problem, since this condition was only introduced 
to prove convergence to reflected Brownian motion, as we discussed above, 
and the local time spent at the first ko rows is negligible. 



3.3 Transition probabilities and scaling limit 

In the previous subsection we described in words how one should define the 
transition probabilities of the uniform walk in the generic wedge W = W a s. 
Here we shall complete the uniform walk construction by providing explicit 
expressions for all of the transition probabilities. We then complete the 
proof of Theorem ^ by computing the directions of reflection for the RBM 
obtained in the scaling limit. 

To get us started, given a vertex x of G we set Xq = x. In the interior of 
the wedge we define the transition probabilities as before. That is, when X* 
is an interior vertex, X^ +1 is to be chosen from the six neighbouring vertices 
with the probabilities a, b and c as in Figure ^ At the apex, that is, if 
X% = 0, the walk moves to either of the vertices on row 1 of the lattice with 
probability l/iV(l). 

If til and ur are both nonnegative, then it only remains to specify the 
transition probabilities from the boundary vertices. However, when either 
riL or tir is negative a little more work needs to be done near the top of the 
wedge, as explained in the previous subsection. In that case, we recall the 
definition of row ko from the previous subsection. For every k = 1, . . . ,ko — 1, 
we set the transition probability from each vertex in row k to each vertex 
in rows k =L 1 equal to (a + c)/N(k). Furthermore, for every vertex in row k 
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that have both a left and a right neighbour, the transition probabilities to 
these two neighbours are equal to b. To the two vertices on row k on the 
sides of W (that have only one neighbour) we assign transition probabilities 
to the single neighbour and to the vertices themselves equal to b. This takes 
care of all the nonzero transition probabilities near the top of the wedge. 

It remains to define the transition probabilities from the boundary ver- 
tices. As we explained in the previous subsection, we are looking for transi- 
tion probabilities that satisfy conditions 1-3 on page ED Below we shall give 
explicit expressions for these transition probabilities from the left boundary 
vertices for arbitrary ni, > 0. This is in fact sufficient to allow us to derive 
the transition probabilities for all possible wedges (i.e. for all combinations 
of til and tir), as we shall explain first. 

Observe that by left-right symmetry we can derive the transition prob- 
abilities from the right boundary vertices for any given ur, if we know the 
corresponding transition probabilities from the left boundary vertices for 
riL = ur. Because the case til = was already covered in Section 12.21 it 
only remains to show that one can obtain the transition probabilities from 
the left boundary vertices for negative nt from those for positive ul- To do 
so, observe from Figure 0] that the left side of a wedge W with given < 
coincides with the right side of a (different) wedge W with tir = —n^ — 1. 
We propose that at the jth vertex on row k of W one can take the proba- 
bility of a step S equal to the probability of the step —S at the jth vertex 
on row k of W , counted from the right side. Here we exploit the up-down 
symmetry inherent in condition 1 on page ED It follows that it is indeed suf- 
ficient to provide the transition probabilities from the left boundary vertices 
for positive only. 

So let rii > be fixed. To specify the transition probabilities we will 
need some notation. We write po[j, I] for the transition probability from the 
jth vertex on a row k to the Ith vertex on the same row. By p±\j, I] we 
denote the transition probability from the jth vertex on a row k to the Ith 
vertex on the row k ± 1 (vertices on a row k are numbered 0,1,..., N(k) — 1 
from left to right). If hr and til are both nonnegative, then these transition 
probabilities are to be used at all rows k > of the lattice, otherwise they 
are valid for the rows k > k$. Finally, we write qj[S] for the probability 
of a step S from the jth vertex on a row k. Remember from the previous 
subsection that there are Nl = ni + 1 left boundary vertices on the rows 
of G, so that we have to give the transition probabilities for j = 0, 1, . . . , til- 

First we specify the transition probabilities from a given row to the row 
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above. For all j = 0, 1, . . . , hl, 

V- [j, 0] = qj [u+(u + v) (n L - j) + ih] = —^—a. (30) 

Second, for all j = 0, 1, . . . , ul the transition probabilities to the same row 
are given by 

Po\j,3 + l]=q j [u + v] = ^-b, (31) 
PotiJ " 1] = <&[-(« + «)] = — XT^ (32) 

ra L + 1 

Pob. J] = ?;P = — — — 1 — 6 - ( 33 ) 

™L + 1 

Third, specifying the transition probabilities to the row below is a bit more 
complicated. For each of the boundary vertices j = 0, 1, . . . , ul we have 

P+ [j, 2n L + 1] = g ■ [-u + (u + v)(n L + l-j)- ih] = . (34) 

TIL + 1 

However, for j = we have 

p+[0,n L J = qj[-u - ih\ = — — , (35) 

nL + 1 

p+[0, 0] = qj[—u — (u + v)nx — i/i] = a + c, (36) 
whereas for j = 1, 2, . . . , ul — 1, 

(a + c)n£ 

P+b> jj = ?7 -« - (« + - in = — j— , (37) 

nL + 1 

P+ [j, «L + j] = qj [-u- ih] = ( Q ± C ^ L j ( 38 ) 

«l + 1 

ci H - c 

P+b'> 2 i] = - (« + f)K - j) - ih] = — , (39) 

n L + 1 

and finally, for j = til we have 

p+[n L , 2ul] = qj[~u - ih] = a + c, (40) 
p+[n L , 2{n L - n) - 1] = ^-[-u - (u + v)(2n + 1) - i/i] = - + - - . (41) 

In the last equation, n is an integer taking values in {0,1,..., ul — 1}- 
The above list specifies all the nonzero transition probabilities from the left 
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Figure 5: Transition probabilities for the reflected random walk in a wedge 
with riL = 1 and ur = 0. The inset shows the lattice dimensions. 



boundary vertices. Figure 03 shows an example of the transition probabilities 
in the case ni = 1. 

We deliberately gave both the transition probabilities and the corre- 
sponding step probabilities at the boundary vertices, to make it easy to 
verify that conditions 1-3 on page 1161 are indeed satisfied. To verify condi- 
tions 1 and 2, one simply has to add up the relevant transition probabilities. 
To check condition 3, one may compute the real and imaginary parts of the 
first step of the random walk from each of the til + 1 boundary vertices 
from the step probabilities (the result is given in the following paragraph). 
As we explained in the previous subsection, the random walk defined above 
will converge in the scaling limit to a reflected Brownian motion. Moreover, 
the walk is uniform by construction and satisfies Lemma |2I Thus, the RBM 
will hit any horizontal cross-section of W with the uniform distribution. 

The proof of Theorem ^wiU therefore be complete if we can show that 
the angles of reflection of the RBM are $l = a and "&r = (3. As in Section^ 
the direction of reflection with respect to the left side of W is given by the 
expected value of the first step from a left boundary vertex. From the step 
probabilities we compute 



E 2 [ReSi] 



1 



c(v — u) + b{u + v) 



n L + l 



(a + c)(2niu + n 2 L (u + v)) 



(42) 



E* [Im Si] 



1 



h(2c + 2 (a + c)n L ) 



(43) 



n L + l 
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where z is any left boundary vertex. Observing that u/h = cot <p and v/h = 
cotip, this gives us the following expression for the reflection angle $l with 
respect to the left side of the wedge: 

cot ip[(a + c)(n\ + 2n£) — b + c] + cot ijj[(a + c)n| — b — c] 
cot(i7x + a) = 



2c + 2(o + c)n L 
Substituting Equations ©-© yields 



(44) 



, „ . [(cot ip + cot ib)nr + cot 09I 2 — 1 , 
cot(# L + a = iy , ^ . 45 

According to Equation (|22j). the result simplifies to 

cot 2 a — 1 . . , . 

cot^ L + Q = = cot2o). 46 

2 cot a 

Thus, the reflection angle of the Brownian motion with respect to the left 
side is simply = a. By symmetry, the angle of reflection with respect to 
the right side is (3. This is consistent with the results obtained in Section El 
and completes the proof of Theorem ^ 



4 Properties of the RBMs 

This section reviews several properties shared by the reflected Brownian 
motions with constant reflection angles, and sheds some light on the RBM Q g 
when a+/3 > ir. The properties discussed below will be used in the following 
section to derive distribution functions for the RBMs. 



4.1 Intertwining relations for the uniform walk 

From the definitions of the uniform walks in Sections [21 and 01 the following 
interesting picture arises. Started with the uniform distribution from the 
vertices in a given row of the lattice, the walk can be seen as a walk from 
row to row on the lattice that remains uniform on each row all the time. 
This can be stated more precisely in the form of an intertwining relation, as 
was noted for the case of symmetric wedges by Dubedat (!]. Here we shall 
describe the generalization of his result to generic wedges. 

First let us explain what is meant by an intertwining relation. Let 
{Pt ■ t > 0) and (.P/ : t > 0) be two Markovian semigroups with dis- 
crete or continuous time parameter and corresponding state spaces (S,S) 
and (S',S'). Suppose that A is a Markov transition kernel from S' to S. 
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That is, A is a function A : S' x S — > [0, 1] such that (1) for each fixed 
x' € S' , A(x' , • ) is a probability measure on S, and (2) for each fixed A E S, 
A( • , A) is <S'-measurable. The two semigroups (Pj) and (P/) are said to be 
intertwined by A if for all t > and every pair (a', A) € 5' x 5, 



In a more compact notation, the semigroups are intertwined if for all t > 
the identity APt = P/A between Markov transition kernels from S' to S 
holds. Examples of such intertwining relations have been studied before, see 
for instance [21 E3 • 

In our case, we are interested in the uniform random walk (X n ) on the 
graph G = G a ^ covering the wedge W = W a ^, as defined in Section |3J 
Let (P n ) be its semigroup, and let (P^) be the semigroup of the random 
walk on the set N of natural numbers with transition probabilities 



for k > and p'(0, 1) = 1. Observe that p'(k,l) is just the conditional 
probability that X n+ \ will be on row / of G, given that X n is uniformly 
distributed on the vertices of row k. 

Now consider the Markov transition kernel A from N to G such that for 
each k € N, A(k, • ) is the uniform measure on the vertices of row k of G. 
Recall that for the walk (X n ), if Xq is uniformly distributed on row k, the 
walk will stay uniform on the rows of G afterwards (Lemma EJ)- It follows 
that for each k G N, A C G and n > 0, 



Hence, we have the intertwining relation AP„ = P^A. The reader should 
compare this with the statement and proof of Lemma |2j 

This discrete intertwining relation may be used to compute the Green 
function for the uniform walk (X®), killed when it reaches the row M > 0. 
The computation of this Green function gives us the expected local time 
spent by the walk on a given boundary vertex of G, and thereby proves 
Lemma |31 

Proof of Lemma OH Consider the Green function G' M = (I — P^) _1 = 
Y^n=o{P'M) n f° r t ne random walk on N with transition matrix P' M as in- 
troduced in Equation (|48j) . except that now the walk is killed as soon as it 




(47) 




(48) 



J2Mk,x)Pn(x,A) = J2K(kJ)Hl,A). 



(49) 
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reaches level M. It is not difficult to verify that the first row of this Green 
function is given by 



{N{1) \ a+c) N(k) forfc>0; 
G' M [0,k} = { (50) 
1 + WV)9mM forfc = 0. 

Now consider the Green function Gm for the uniform walk (X®) on the 
graph G in W, killed when it reaches {z : Imz = —Mh}. Let us denote the 
jth vertex on row k of G by (j, k). Then, by the intertwining relation (|49|) 
(or by Lemma |2J), the Green function Gm is related to G' M by 

G M [(0,0),(j,fc)] =^—G' M [0,k) for j = 0,l,...,iV(fc)-l, (51) 

since the expected local time spent at vertex (j, k) by the walk (X®) before 
it is killed at row M is the same for all j = 0, 1, ... , N(k) — 1. It follows 
that the expected local time spent at any vertex (j, k) by the walk before it 
is killed at row M is smaller than (a + c) . □ 



4.2 Intertwining relations and time reversal of the RBMs 

In the previous subsection we considered the intertwining relation between 
the uniform walk on G and a random walk on the integers. Here we will 
turn our attention to the scaling limit. This time, let (Zt) be the RBM Q ( g in 
W = Wq.,/3, and let (Pt) be its semigroup. Consider the Markov transition 
kernel A from the positive reals R+ to W which for each fixed y £ R+ assigns 
the uniform measure to the horizontal interval [— y cot a — iy,y cot (3 — \y\. It 
is clear from the intertwining relation (|49fl between the random walks, that 
(Pt) and the semigroup of the scaling limit of the random walk on N will be 
intertwined by A. It remains to identify this scaling limit. We claim that this 
is a 3-dimensional Bessel process, or in other words, it is a Brownian motion 
on 1R, + conditioned not to hit the origin. This generalizes Proposition 1 in 
Dubedat |H to the following statement: 

Theorem 5 Let (Z t ), (Pt) and A be as above and let (P() be the semigroup 
of the 3-dimensional Bessel process taking values in Then (Pt) and (P/) 
are intertwined by A. In particular, for all y > 0, the process (^ — ImZ^ y ' ^ 
is a 3-dimensional Bessel process started from y. 
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Proof: As we remarked above, it is sufficient to identify the scaling limit 
of the random walk {X' n ) on N. Since the rows of the lattice have spacing h, 
the proper scaling limit is obtained by considering the linear interpolations 
of the processes hX', 2t Jna where a 2 = 2(a + c)h? is as before, and then 
taking n — > oo. We may derive the infinitesimal generator for the limit by 
computing, for a sufficiently differentiable function / on R+, 

" + ' ( f (x + —) N(nax/h + l)-f( x -—) N(nax/h - 1] 



N(nax/h) \ \ na J \ na 

+ 2bf(x) - f(x) = - 2 (-f'(x) + lf"(x)) + o(n- 2 ). (52) 
n z \x 2 J 

Here, 1/n 2 is the time scaling, and we recognize the generator of the 3- 
dimensional Bessel process (see Revuz and Yor |12| Chapter VI, §3 and 
Chapter III, Exercise (1.15) for background on the 3-dimensional Bessel 
process and its semigroup). □ 

We now turn our attention to the time-reversal of the RBMs. Precisely, 
let (Zt) be an RBM a a in the triangle T = T a n started from the top, and 
stopped when it hits [0,1]. We are interested in the time-reversal (Zt) of 
this process. From the time-reversal properties of the 3-dimensional Bessel 
process [121 Proposition VII(4.8)] we know that until (Zt) first hits the 
boundary of T, it is a complex Brownian motion started with the uniform 
distribution from [0, 1] and conditioned not to return to [0, 1]. In fact, we can 
describe the full process (Zt) in terms of a conditioned reflected Brownian 
motion. This is again a generalization of an earlier result of Dubedat [U 
Proposition 2]: 

Theorem 6 The time-reversal of the RBM a ^ in the triangle T = T a p, 
started from the top and stopped when it hits [0, 1], is the RBM n _ a n _p in T 
started with the uniform distribution from [0,1], conditioned not to return 
to [0, 1], and killed when it hits the top of the triangle. 

Proof: Recall the Green function of the uniform walk (X®) of Equa- 
tion (|51j) . By Nagasawa's formula (Rogers and Williams |14| III. 42]) the 
time-reversal of this random walk is a Markov process with transition prob- 
abilities 

u- u\ n m Gj\/[(0,0), (l,m)]p[(l,m), (j,k)] { 

qM[(j,k),(l,m)\ = , . 53 

Gm (0, , (j, k)\ 
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Here, the p[(l , m) , (j , k)] are the transition probabilities for the walk (X®) 
killed at row M, as specified in Section 13.31 We are interested in the tran- 
sition probabilities for the reversed process in the limit M — > oo. From the 
expression for the Green function it is clear that in the limit one gets for 
k, m > 

q[(j,k),(l,m)} =p[(l,m),(j,k)}. (54) 

Observe in particular that in the interior of the wedge we recover the tran- 
sition probabilities of the original walk. Hence, the reversed walk converges 
to a Brownian motion in the interior. 

Moreover, by condition 1 on page 1161 it is clear that at every vertex on 
any row k > of the lattice, the probability that the reversed walk will 
step to the row k + 1 is equal to the probability to step to the row k — 1. 
Note especially that this is not only true at the interior vertices but also 
at the boundary vertices. Therefore, the expected step of the random walk 
from any vertex on the rows k > is real. In particular, on the sides of W 
the walk must receive an average reflection in the real direction toward the 
interior of W. It follows that the time-reversal of the reflected Brownian 
motion in the wedge is a reflected Brownian motion with reflection angles 
tt — a and tt — j3 with respect to the sides of W. □ 

4.3 Reflected Brownian motions for a + (3 > n 

The fact that the RBM Qj/ g is intertwined with a three-dimensional Bessel 
process sheds some light on the behaviour of the reflected Brownian motions 
with reflection angles satisfying a+fi > tt. It is the purpose of this subsection 
to look at these RBMs more closely. For the duration of this subsection we 
will fix a, [3 £ (0, tt) such that a + j3 > tt. Then an RBM Qij a in the domain 
T = T a> p may be described by considering an RBM 7r _ Q 7r _ p in the wedge 
W^-a^-p and putting it upside-down, as we shall explain below. 

We set x := cos a sin/3/ sin(a + /3) and z := — sin a sin f}/ sin(a + (3). Let 
(Z t ) denote RBMn-a^-p in the wedge W^-a^-p, and consider in particular 
the process (^Z^ v+Z ' with y € R+. Here, A is the Markov transition 
kernel from R+ to W^-a^-p introduced in the previous subsection. Let / 
be the transformation f : w *— > w + x — iz which puts the wedge upside-down 
as illustrated in Figure El and set Y t y := f [z^^'^ . Then the process 

(Xt)i stopped when it hits the interval [0, 1], is an RBM a/ g in T a ^ started 
with the uniform distribution from the horizontal line segment at altitude y. 

From the intertwining relation of Theorem we conclude that the imag- 
inary part of (Y^) is a Brownian motion conditioned not to hit —z. It 
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Figure 6: By putting a wedge upside-down we can shed some light on the 
RBM Qij a when a + (3 > ir. 

follows (see Corollary VI(3.4)]) that (Y^) has positive probability to 
never reach [0,1]. In fact, the probability that (Y^) does reach [0,1] is 
z/(z + y). It is furthermore clear from Theorem [5] that given the event that 
(Yt) does reach [0, 1], it will arrive there with the uniform distribution. 

Now let P y denote probability with respect to the process (Y^). Then 
by what we said above, (y/z + 1) P y is a probability measure on reflected 
Brownian paths in T Qt g started from the horizontal line segment at altitude y 
that end on [0, 1]. Taking the limit i/->oowe obtain a conformally invariant 
probability measure on paths of RBMq, r in T a a that start "with the uniform 
distribution from infinity" and arrive on [0, 1] with the uniform distribution. 
Henceforth, when we speak about RBM ajj g with a + (3 > tt, we shall always 
assume that we restrict ourselves to this collection of Brownian paths and 
the corresponding probability measure introduced above. 

4.4 Conformal invariance and locality 

Two elementary properties shared by the RBMs are conformal invariance 
and the locality property. To explain what we mean by these properties, 
let us first consider an RBM Q)( g in the wedge W = W a ^ started from the 
point x € W. Call this process {Zf ), and let v\ := exp((2a — 7r)i) and 
V2 := exp(— 2(3 i) denote the reflection fields on the left and right sides of W. 
Then it is well-known (compare with Equation (2.4) in |15j ) that we can 
uniquely write 

Zf = Bf + v 1 Y t l +v 2 Y t \ (55) 

where (Bf) is a complex Brownian motion started from x, and (l^ 1 ) and (Y^) 
are real-valued continuous increasing processes adapted to {Bf ) such that 
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Figure 7: The locality property says that the RBM aj/ g in T started from w 
and stopped when it exits from the subset D behaves just like an RBM a ^ 
in D started from w. 



Yq = Yq = 0. Moreover, Y^ increases only when Zf is on the left side of W 
and Yj 2 increases only when Zf is on the right side of W (Y^ and Y 2 are 
essentially the local times of {Zf) on the two sides of W). 

Now let g be a conformal transformation from W onto a domain D with 
smooth boundary. Consider the sum 



N-l 



g(Zf) - g(ZS) = £ [ 9 (Zf j+1)t/N ) - g(Zf t/N )] . (56) 

3=0 

To compute the sum we use Taylor's theorem to expand each term. The 
computation is very similar to the one in Section [2.31 In particular, on the 
boundary we may keep only the first-order terms. Then letting N — > oo and 
using the fact that the real and imaginary parts of g are harmonic, just as 
in the proof of Ito's formula (see e.g. Sections 4.2 and 4.3 in Gardiner [Jj for 
a nice discussion) one obtains 

g(Z*)-g(ZS)= [ t g'(Z-)dB-+v 1 f g\Z*) dY s 1 +v 2 f g> \Zf) dY s 2 . (57) 
Jo Jo Jo 

The first integral in (|57|) is the usual expression for the conformal image 
of Brownian motion. Making the usual time-change u(t) := Jq \g' (Zf)\ 2 ds 
(see Revuz and Yor ^1 Theorem V(2.5)]) and denoting its inverse by t(u), 
we conclude from Equation (fSTj) that the process Z u := g(Zf^) is a re- 
flected Brownian motion in D with reflection vector fields v\ g'{g~ l { ■ )) and 
V2 g'(g~ 1 ( ■ )) on the two "sides" of D (i.e. the images of the two sides of W). 
Note in particular that because g is an angle-preserving transformation, the 
process (Z u ) is also reflected at the angles a and (3 with respect to the 
boundary of D. This shows that the RBM Qj/ g is conformally invariant. 
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We may use the same reasoning to explain what we mean by the locality 
property of the RBM Q «. We change the setup to one that will be more 
useful later. Indeed, we now let (Zt) be an RBM Qj/ 3 in the triangle T = T Qj/ g 
started from the top w = w a ^, and set r := inf{i > : Z t £ [0,1]}- 
Furthermore, we take g to be a conformal map of T onto an open subset D 
of T that fixes w, and such that the left and right sides of T are mapped 
onto subsets of themselves. Finally, we set a := inf{t > : Zt € T \ D}, the 
exit time of the RBMq, ^ from the subset D. See Figure[7|for an illustration. 

From the calculation above we conclude that up to the stopping time r, 
the process g(Zt) is just a time-changed RBM a ( g in the subset D of the 
triangle T. In particular, modulo a time-change the laws of (g(Zt) : t < t) 
and {Z t : t < a) are the same. This is what is called the locality property. 
For more background and for consequences of the locality property we refer 
to the SLE literature jHllTfi]. 

5 Distribution functions 

In this section we compute several distribution functions associated with the 
family RBM ai/ 3 of reflected Brownian motions. We fix a and (3 in (0, tt) for 
the duration of the section, with no further restrictions on a, (3 (recall that 
when a + (3 > tt, we assume that we work with the probability measure 
of Section f4..'-{|h Furthermore, we fix two angles A, fi € (0,7r) such that 
A + fi < tt. These two angles define the domain T = T\ ^ in this section. 
To be somewhat more general, we will derive distribution functions for the 
RBM a $ in the triangle T = T\„ rather than in the domain T a ^. 

5.1 Notations 

First we introduce some notations. We shall call the set of points discon- 
nected from {0, 1} in the triangle T by the path of the RBM up to the 
time when it hits [0, 1] the hull K of the process. This hull has exactly one 
point in common with the real line, which we denote by X. We shall denote 
by \w\Y the distance of the lowest point of the hull on the left side to the 
top w = wa,/ii and by \w — 1\Z the distance of the lowest point of the hull 
on the right side to the top w. Thus, all three random variables X, Y and Z 
take on values in the range [0,1]. See Figure |H] for an illustration of the 
definitions. 

Below we shall compute the marginal and joint distributions of the vari- 
ables X, Y and Z. We shall see that these can be expressed in terms of the 
conformal transformations of the upper half of the complex plane onto the 
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Figure 8: Definition of the hull K of the reflected Brownian motion with 
parameters a and (3 in the triangle T, and of the random variables X, Y 
and Z. 



triangles T 7j( 5. Thus it is useful to review some properties of these trans- 
formations first. We simplify the notation somewhat by writing 7' as an 
abbreviation for 7/V whenever 7 denotes an angle. 

By the Schwarz-Christoffel formula of complex analysis (see ^ Sec- 
tion 6.2.4] or Section XI.3]), the unique conformal transformation of 
the upper half-plane H onto that fixes and 1 and maps 00 to w 1) s is 
given by 

where B(j',5') = F(Y)T(S')/T(j' + 5') is the beta-function, and B z (y,5') 
the incomplete beta-function (see e.g. |S] for background on these special 
functions). Symmetry considerations show that the transformations satisfy 

F 7 ,i(«) = 1-^(1 -u) and F-j(x) = l-F^(l-x). (59) 

See Figure |H1 for an illustration. 

A different kind of distribution that we can compute is the conditional 
probability that the last side of the triangle visited by the RBM, given that 
it lands at X = x, is the right side. As we shall see, this distribution can 
also be expressed in term of triangle mappings. In fact, it turns out that 
there is a remarkable resemblance between this conditional probability and 
the marginal distribution functions of the variables X, Y and Z. 

5.2 Characteristics of the hull 

Here we derive the (joint) distribution functions of X, Y and Z, which are 
characteristics of the hull generated by the RBM. Remember that we are 
considering an RBM Qij g in the triangle T = T\ tfl started from w = wx^. For 
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Ox 1 u 1-u 1 l-x 1 

Figure 9: Transformations of the upper half-plane onto triangles. 

convenience, let us also introduce the angle v := tt — A — [i. Then our main 
conclusion may be formulated as follows: 

Proposition 7 Let a = a(y) and b = b{z) be the points on the left and right 
sides ofT at distances \w\y and \w — \\z from w, respectively. Then the joint 
distribution of X, Y and Z is given by 

P[X < x,Y < y,Z < z] = (60) 
F / ^>)-^>) \ F ( ~Fx» \ 



where the images of a and b under the map F^ can be expressed in terms 
of y and z as 

1 

F *> ) = i-^(i-*) = W (62) 



^» = i-i^t7t; ( 61 ) 



Note that in the last equation we used the symmetry property (|59|) . 

Proof: The idea of the computation of P[X < x, Y < y, Z < z] is illus- 
trated in Figure 1101 Consider an RBM Qij a in the triangle T started from 
the top w, and stopped as soon as it hits the counter-clockwise arc from a 
to b on the boundary (the thick line in the figure). Then the probability 
P[X < x,Y < y, Z < z] is just the probability that this process is stopped 
in the interval (0, x). 

We now use conformal invariance and locality. Let g = g a (y),b(z) De the 
conformal map of T onto T a ^ that sends a to 0, 6 to 1 and w to u> Qj/ 3, as 
illustrated in Figure ITU1 Then the probability that we are trying to compute 
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Figure 10: This figure illustrates how the joint distribution function of the 
random variables X, Y and Z can be computed. As explained in the text, 
the joint probability P[X < x, Y < y, Z < z] is just g(x) — g(0). 



is exactly the probability that an RBM Q / g in T a ^ started from w a ^ and 
stopped when it hits [0, 1], is stopped in the interval (g(0), g(x)). But since 
the exit distribution of the RBM is uniform in T Qj/ g, this probability is simply 
g(x) - g{0). Thus, 

P[X <x,Y <y,Z <z]=g(x)-g(0). (63) 

Our next task is to find an explicit expression for this joint probability by 
deriving the explicit form of the map g = g a i y )u z )- The explicit form of g is 
obtained by suitably combining conformal self-maps of the upper half-plane 
with triangle mappings. How this is done exactly is described in Figure ITT1 
The expression for the joint distribution follows. □ 



By sending one or or two of the three variables x, y and z to 1, and 
using the symmetry property (|59[). one may derive the following corollaries 
of Proposition 

Corollary 8 We have the following joint distributions: 

P[X<x,Z<z]= F a ,p (F A -J(x)F-i(*)) ; (64) 
P[X < x ,Y<y]= F P;a - F p>a (F-{(1 - x)F£(yj) ; (65) 



P[Y < y, Z < z] = F a>f3 ( Jp\ [z) , V^TT-A (66) 

V Fu,,m (z) + F u X (y) - F V4l (z)F^ x (y) J 

+ F ( F - 1{V) \ 1 

P ' a \F^{z) + F-\ (y)- F^ (z)F- l x (y)J 
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Figure 11: This illustration shows schematically how one obtains an explicit 
form for the map g in terms of the variables y and z. The notations a, b 
and x in the figure are short for F^(a), F^(b) and FT^(x). 



Corollary 9 The marginal distributions of X , Y and Z are given by 

P[X<x] = F a ^(F^{x))- (67) 

P[Y<y] = Fp, a (F-{{y))- (68) 

P[Z<z] = F a ^(F~l{z)). (69) 



Observe that the marginal distributions of the variables X, Y and Z take 
on particularly simple forms. These marginal distribution functions have a 
nice geometric interpretation. For instance, P[Y < y] is just the image 
of y under the transformation that maps the triangle T u ^\ onto Tp jCt , fixes 
and 1 and takes w v ^\ onto wp^ a . Similar observations hold for P[X < x] 
and P[Z < z). 

These observations lead to some intriguing conclusions. First of all, we 
conclude that for an RBM^ ^ in an equilateral triangle (so that a = (3 = 
A = n = v = 7r/3) all three variables X, Y and Z are uniform. This is 
not so surprising when we realize that the hull in this case is the same as 
that of the exploration process of critical percolation, as we noted in the 
introduction. Indeed, in the case of percolation X and Y can be interpreted 
as the endpoints of the highest crossing of a given colour between the sides 
(0, w) and (0, 1) of the triangle. Thus by symmetry, if X is uniform, then 
so are Y and Z. 
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Figure 12: The unions Kq and K' of the shaded sets on the left and right 
are generated by different RBMs, but have the same law. 

In other triangles, similar but more intricate connections exist. For 
example, let (Z t ) be an RBM aj/ g in T = T\ ijU started from w = w\ tfM 
and stopped when it hits [0, 1] , as before. Compare this process with an 
RBM (Z' t ) in T started from 1, stopped on [0,w] and reflected on (0,1) at 
an angle a and on (w, 1) at an angle with respect to the boundary. Here it 
is assumed that small angles denote reflection away from 1. To this second 
RBM we can associate normalized random variables X', Y' and Z', measur- 
ing the distances of the exit point on [0, w] to w and of the "lowest points" 
of the hull on [w, 1] and [0, 1] to w and 0, respectively. See Figure It 
follows from Corollary El that X and Z' have the same distribution (and so 
do Y and X'). 

This result has an interesting interpretation in terms of the hulls gen- 
erated by the two processes, as we shall now describe. We write VL for the 
collection of closed, connected subsets C of T such that the right side of T 
is in C and T \ C is connected. We further define Q as the collection of 
compact A C T such that A = A Pi T, and T \ A is simply connected and 
contains the right side of T. We then endow with the <7-field generated by 
the events {K € SI : K n A = 0} for all A € Q. This setup is similar to the 
one in Lawler, Schramm and Werner [01 Sections 2 and 3]. In particular, a 
probability measure P on O is characterized by the values of P [K n A = 0] 
for A G Q, see Lemma 3.2]. 

Theorem 10 Consider the processes (Zt) and (Z' t ) stopped on [0, 1] and 
[0,w] as described above. Let Kq and K denote the sets of points in T that 
are disconnected from by the paths of (Zt) and (Z' t ), respectively. Then the 
laws of Kq and K on the space fi are the same. 
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Proof: Let P be the law of Kq, and let A G Q. Denote by a and b the 
points of A n dT closest to w and 1, respectively. Let g : T —* T \ ^4 be the 
conformal transformation that fixes 1 and w, and maps onto a if Ima > 0, 
and onto otherwise. Then, by conformal invariance of the RBM a a in T, 
P[K n A = 0] = PLY G (c/- 1 (6), 1)] (compare this with our discussion of 
the locality property in Section I4.4JI . Likewise, the law P' of K'q satisfies 
P'[K n A = 0] = P'[Z' G (^(fr),!.)]. Since X and Z' have the same 
distribution by Corollary the theorem follows. □ 

5.3 Last-visit distribution 

As we discussed in Section 14. 2| the imaginary part of the KBM a ^ in T Qj/ g 
stopped when it hits [0, 1] is a three-dimensional Bessel process, and so is 
the imaginary part of its time-reversal. In particular, the time-reversed pro- 
cess, considered up to the first time when it hits the left or right side of the 
triangle, is a complex Brownian motion started with uniform distribution 
from [0, 1] and conditioned not to return to the real line. In other words, up 
to its first contact with the left or right side, this process is a Brownian ex- 
cursion of the upper half-plane, started with uniform distribution from [0, 1] 
(for background on Brownian excursions, see [B] and [H|)- This fact allows 
us to derive the following result: 

Proposition 11 Let (Zt) be an RBM a p in the triangle T = T\„ started 
from w = wx,^. Let r := inf{t > : Zt G [0, 1]}, and let a be the last time 
before r when (Zt) visited the boundary of T . Let E denote the event that 
Z a is on the right side of the triangle. Then 

P[E | Z T = x] = i^-a,*-/? ■ 

It is shown in Dubedat jlj how one derives this result in the special case 
where a = f3 = \ = fi = ir/3. Here, for the sake of completeness, we repeat 
the computation for the general case. 

Proof: We want to use the fact that the time-reversal of the RBM Qj/ g in 
the triangle T ai p starts out as a Brownian excursion. So we first map T\ :fl 
onto T Qj /3 by the transformation F a ^ o , This maps the point x onto 

y-=F a , p (f^(x)). 

Next, let (Yt : t > 0) be a Brownian excursion of the upper half-plane, 
and let S be the first time when (Yt) visits either the left or right side of T a p. 
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Figure 13: A complex Brownian motion Bt started from y + ie in the tri- 
angle T a; p. This process conditioned to exit the strip {z : < Imz < M} 
through the top boundary and stopped when it hits the boundary of the 
triangle T a n corresponds in the limit M — ► oo, e { to the time-reversal of 
the RBM Qj/ 3 in the triangle, as explained in the text. 



Then 

P[E | Z T = x) = lim / P y+is [Y s e dz], (70) 

where the integrals runs over the right side of the triangle T a p, and P z 
denotes probability with respect to the Brownian excursion started from z. 

Now let (Bt : t > 0) be a complex Brownian motion, let U be the first 
time when (Bt) visits either the left or right side of T a p, and let Um be 
the first time when (Bt) exits the strip {z : < Imz < M}. Suppose that 
P' z denotes probability with respect to the Brownian motion started from z. 
Then, using the strong Markov property of Brownian motion (see Figure ITUl 
for a sketch), we have 

P y+i£ [Ys G dz] = lim P' +i£ [B u edz\ImB UM =M] 

M — >oc 

= P'JBu e dz lim , zl r ^ 

V+1£l l M^P' y+i£ [lmB UM =M] 

Im z 

= p; +iE [%Gdz]— , (71) 

where in the last step we have used j!21 Proposition 11(3.8)]. Combin- 
ing Equations (|70)) and l|71jl. we see that we have to compute the limit of 
P' y+ie [Bu € dz]lmz/e as e — > 0. This computation can be done by using 
the conformal invariance of Brownian motion. 

Remember that the probability that a complex Brownian motion started 
from a + \b leaves the upper half-plane through (— oo,x), is given by the 
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Figure 14: Conformal transformations between the upper half-plane and the 
domains T a ^ and T^-a^-p- 



harmonic measure u{x) of (— oo,x) at the point a + ib. It is straightforward 
to verify that 

f x b dt 11 x — a 

LO(x) = / —-n ; ttt = — I arctan — - — . 72 

v ; i_oo vr b 2 + (t - a) 2 2 7T b y 1 

Thus, mapping the triangle T a ^ conformally to the upper half-plane by the 
transformation F~ n as in Figure ll~H we can write 

UmP^edz]— = wKAa)(F - i{z) _ a)2 P3) 

where a = F~^(y) = (x). Therefore, using (|58|). 

P[£7 | Z T = x] (74) 

1 Z" 00 dt T „ . . 

■lmF ai p(t) 



^-"'(l-o) 1 -^' f 00 dt 



7T Ji (t — ay Ji 



sin/3 !_ a / 
a 


(1- 


-a] 


I 1 " 


-ff 


L 


7T 












sin^ x _ / 
a 


(1- 


-a] 


I 1 " 


-0' 


L 


7T 












sin^ , 




-a, 


I 1 " 






a 


(1- 




L 


7T 













dt 



u (t- a) 



2 



3C 



« a '- 1 («-l)^- 1 (u-a)- 1 d« 

i 



where in the last step we have made the substitution t = l/u. 

Using equations 15.3.1 and 15.2.5 for the hyper geometric function in jllj 
and the formulas 6.1.15 and 6.1.17 for the gamma function from jS], we 
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finally derive 
P[E \Z T = x] 



(75) 



smf3T(2-a'-P')T(p') x _ a , 




-P' 2 F 1 (l,2-a' -/3';2-a';a) 





= F 7 




This is the desired result. 



□ 



In words, we have considered the conditional probability that the last 
side visited by an RBM a ^ in started from w\ jfM is the right side, given 
that the exit point X equals x £ (0,1). This conditional probability is 
exactly given by the image of x under the transformation that maps T\ „ 
onto T^-u^-p, fixing and 1 and sending onto vj^ q,7t /3 • 
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